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ABSTRACT  OF  THE  THESIS 


A  Statistical  Survey  of  Electron  Pitch  Angle  E>istributions 
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by 

Timothy  F.  Alsruhe 
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S.  V.  Venkateswaran,  Chair 


The  diffuse  aurora  is  produced  by  precipitating  electrons  whose  origin  is  believed  to  be  the 
central  plasma  sheet.  The  process  by  which  these  electrons  are  continually  scattered  into 
the  atmospheric  loss  cone  is 'apparent!^  unresolved.  Arguments  included  processes  which 
cause  pitch  angle  scattering  at  the  strong  diffusion  limit.  A  survey  of  equatorial  pitch  angle 
distributions,  of -rO.2-20KeV  electrons  has  been  conducted  to  determine  the  extent  of  pitch 


angle  diffusion.  Pitch  angle  distributions  have  been  binned  by  location  and  geomagnetic 
activity  for  the  nightside  near  earth  plasma  sheet.  The  average  pitch  angle  anisotropy  has"^ 
^^^bccn  determined  for  each  bin.  Representative  samples  were  selected  from  the-bins-of" 
geomagnetically  quiet  (Kp<3)  distributions.  The  rate  of  diffusion  was  estimated  by  fitting 

y  .  ^ 

distributions  to  plots  of  theoretically  diffused  distributions.  Of  the  three  energies  examined 
(1.09, 4.52, 19.4keV),  only  distributions  consistent  with  moderate  to  weak  diffusion  were 
observed.  Even  test  cases  oTentirely  isotropic  distributions^ere  consistent  with  only 
moderate  diffusion.  These  results  lead-mer^  conclude  that  strong  diffusion  is  usually  not 
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in  effect  in  the  quiet  near  earth  plasma  sheet.  Electron  Cyclotron  Harmonic  (ECH)  waves 

^  i..  . 

at  the  3/2  harmonic  have^^»e^n  suggested  as  the  scattering  mechanism  responsible  for 
maintaining  the  diffuse  aurora.  The  necessary  ECH  wave  amplitudes  required  to  recreate 
average  distributions  were  calculated.  Derived  amplitudes  were  one  to  three  orders  of 
magnitude  greater  than  those  amplitudes  observed.  Thus,  it  is  concluded  ECH  do  not 
occur  with  the  frequency  and  amplitude  necessary  to  sustain  the  continuous  diffuse  aurora. 
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INTRODUCTION 


The  magnetosphere  as  a  conduit  for  solar  particles  entering  the  earth's  ionosphere  is 
well  established.  It  is  generally  agreed  the  source  of  the  ionosphere's  diffuse  aurora  is  the 
near  earth  central  plasma  sheet  [Lui  et  al.,  1977;  Meng  et  al.,  1979].  The  process  by  which 
electrons  are  scattered  into  the  atmospheric  loss  cone  and  precipitated  in  the  diffuse  aurora 
is  still  a  source  of  controversy.  Kennel  et  al.  [1970]  first  suggested  electron  cyclotron 
harmonic  (ECH)  waves  as  a  mechanism  for  scattering  electrons  at  the  strong  pitch  angle 
diffusion  limit.  Lyons  [1974]  derived  the  ECH  wave  emissions  necessary  to  put  plasma 
sheet  electrons  on  strong  diffusion.  Subsequent  wave  observations  did  not  find  adequate 
occurrences  of  ECH  waves  with  sufficient  amplitude  to  maintain  the  continuous  diffuse 
aurora  [Belmont  et  al.,  1983;  Boeder  and  Koons,  1989]. 

In  this  thesis,  I  present  a  survey  of  electron  pitch  angle  distributions  in  the  nightside 
central  plasma  sheet.  The  purpose  of  this  survey  is  to  determine  the  extent  of  pitch  angle 
diffusion  and  the  necessary  ECH  wave  emissions  required  to  produce  observed 
distributions.  Data  points  are  sorted  into  various  binning  arrangements  based  on 
geographical  location  and  geomagnetic  activity.  Statistically  significant  electron  pitch  angle 
anisotropies  from  each  bin  are  selected  and  used  to  calculate  the  diffusion  and  required 
wave  emissions. 

Chapter  2  contains  a  description  of  the  diffuse  aurora  and  inner  edge  of  the  plasma 
sheet  as  well  as  a  chronology  of  relevant  theories  and  observations.  In  Chapter  3,  the 
theories  and  processes  of  pitch  angle  diffusion  and  wave  particle  interactions  are  outlined. 
The  particle  detector  and  spacecraft  which  provided  the  data  used  in  this  research  are 
describe  in  Chapter  4.  Data  selection,  analysis,  and  observations,  including  sources  of 
ambiguity  and  statistical  uncertainties,  are  discussed  in  Chapter  5.  Conclusions  and  areas 
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for  future  research  are  contained  in  Chapter  6. 
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CHAPTERZ 

BACKGROUND 


2.1  THE  DIFFUSE  AURORA 

The  diffuse  aurora  is  a  broad  band  of  steady  auroral  glow  that  exists  continually  on 
the  auroral  oval.  The  diffuse  aurora  illumination  is  produced  primarily  by  precipitating 
electrons.  Mende  and  Eather  [1976]  estimate  that  less  than  30%  of  the  light  is  produced  by 
protons.  The  diffuse  aurora  is  a  constant  feamre  and  is  present  at  all  levels  of  geomagnetic 
activity.  It  is  the  bright  arc  structure  of  the  discrete  aurora  that  the  casual  observer  refers  to 
as  the  "aurora".  Unlike  the  discrete  aurora,  the  diffuse  aurora  is  not  readily  detected  by  the 
human  eye.  Its  true  immensity  is  only  revealed  by  sensitive  spacebomc  photometric 
devices  and  charged  particle  detectors.  Despite  its  mundane  appearance,  the  diffuse  aurora 
is  the  powerhouse  of  auroral  phenomena,  generating  80%  of  the  total  energy  input  into  the 
polar  region  at  solar  maximum  and  no  less  then  50%  at  solar  minimum  [Whalen,  1985]. 

This  paper  concentrates  exclusively  on  electron  precipitation  on  the  nightside.  The 
general  spectral  range  for  electrons  that  are  the  primary  cause  of  the  diffuse  aurora  lies 
between  a  few  hundred  electron  volts  (eV)  and  ~10keV.  Hardy  et  al.  [  1985]  computed  the 
total  hemispheric  electron  energy  flux  as  varying  from  9.98  x  lO^^ergs/s  ster  at  Kp=l  to 
2.06  X  10*^ergs/s  ster  for  Kp>6-.  The  total  hemispheric  electron  number  flux  varies  from 
9.4  X  lO^el/s  ster  at  Kp=0  to  6.41  x  10^*el/s  ster  for  Kp  ^6-.  This  number  flux  can  be 
misleading  since  73-90%  of  the  energy  flux  is  carried  by  only  a  small  number  of  electrons 
whose  energies  are  greater  than  660eV.  The  number  flux  of  these  electrons  is  1.26  x 
lO^^el/s  ster  at  Kp=0  to  only  1.9  x  lO^^el/s  ster  at  Kp>6-..  Contributions  by  particles 
greater  than  lOkeV  decrease  smoothly  [Meng  et  al.  1979;  Kennel  and  Ashour-Abdalla, 
1982;  Lyons  and  Fennell,  1986].  Discussion  in  the  literature  of  the  pitch  angle  diffusion  in 
the  magnetic  quiescent  diffuse  aurora  is  usually  in  the  context  of  describing  other  more 
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dynamic  auroral  phenomena  [Winningham,  1975;  McDiarmid  and  Budzinsld,  1968;  Qiase, 
1968].  What  can  be  gleaned  from  these  observations  is  that  the  electron  precipitation 
creating  the  diffuse  aurora  is  usually  isotropic  or  nearly  isotropic  in  pitch  angle.  This 
isotropy  has  lead  investigators  to  assume  a  filled  atmospheric  loss  cone  in  the  plasma  sheet. 
Kennel  [1969]  notes  a  filled  loss  cone  as  an  indication  of  strong  diffusion  (discussed  in 
detail  in  Chapter  3). 

On  the  nightside,  the  diffuse  aurora  extends  from  dusk  to  dawn  in  a  band  several 
hundred  kilometers  wide.  The  poleward  edge  of  this  aurora  is  weakly  defined.  Lui  and 
Anger  [1973]  determined  a  globally  average  boundary  location  of  68°.  The  equatorial  edge 
of  the  diffuse  aurora  varies  as  a  function  of  geomagnetic  activity  [Gussenhoven  et  al., 
1981,1983  and  references  therein].  The  equatorward  boundary  is  highly  correlated  to 
variations  in  the  Kp  index.  This  is  particularly  true  on  the  nightside.  Gussenhoven  et  al. 
[1981]  found  the  minimum  nightside  equatorward  boundary  to  be  ~71®  at  18()0MLT  in 
quiet  times  (Kp  =  0).  The  boundary's  maximum  extent  during  geomagnetically  active 
times  reaches  57°  at  0300  MLT  for  Kp  =  5. 

Since  the  flux  of  electrons  into  the  diffuse  aurora  is  approximately  twice  that  of 
ions,  the  charge  imbalance  would  be  expected  to  create  parallel  electric  currents.  But  any 
incipient  parallel  current  is  neutralized  by  a  combination  of  cold  ionospheric  electrons 
scattered  onto  magnetic  field  lines  by  incoming  magnetospheric  electrons  and  by 
backscattered  magnetospheric  electrons.  Energies  of  secondary  and  backscattered  electrons 
range  between  a  hundred  eV  and  a  few  keV.  When  parallel  electric  currents  do  occur  on 
field  lines  they  accelerate  electrons  to  create  discrete  structure  in  the  aurora.  Evidence  of 
acceleration  is  a  peak  in  the  electron  distribution.  Though  both  the  discrete  and  diffuse 
auroral  particles  share  a  common  source  --  the  plasma  sheet,  the  diffuse  aurora  is  defined  as 
the  precipitation  of  unaccelerated  particles  and  no  peak  in  the  distribution  is  observed. 

Similarities  in  the  diffuse  aurora  and  the  plasma  sheet  energy  spectrum  lead 
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numerous  researchers  to  conclude  a  connection  between  the  two  regions  [Winningham  et 
al.,  1975,  Eather  et  al.,  1976,  Lui,  1977,  Meng  et  al.,  1979,  Slater  et  al.,  1980;  Howitz  et 
al.,  1982].  Meng  et  al.  [1979],  the  most  widely  referenced  study  on  this  connection, 
examined  observations  from  the  low  orbiting  DMSP  satellite  just  above  the  diffuse  aurora 
and  ATS  6  satellite  in  the  equatorial  plasma  sheet.  Simultaneous  observations  taken  by 
these  satellites  when  on  conjugate  field  lines  produced  nearly  identical  electron  spectra. 

2.2  THE  PLASMA  SHEET 

The  boundaries  of  the  inner  edge  of  the  plasma  sheet  are  more  easily  observed  than 
explained.  Gussenhoven  et  al.  [1981]  projected  the  equatorward  boundary  of  the  auroral 
oval  back  onto  the  magnetic  equatorial  plane.  They  found  the  latitudial  extent  of  the  inner 
edge  varies  in  local  time  and  is  approximately  related  to  the  Kp  index.  Mapping  these 
boundaries  to  field  lines  in  the  equatorial  plane,  the  inner  edge  ranged  from  a  minimum  of 
3.5  earth  radii  (Re)  at  0300MLT  for  Kp  =  5  to  a  maximum  of  10.2Rc  at  21()0MLT  for  Kp 
=  0. 

Two  theories  predominate  on  what  process  is  responsible  for  defining  the  inner 
edge  of  the  plasma  sheet.  One  theory  is  based  on  pitch  angle  diffusion  [Kennel,  1969; 
Ashour-Abdalla  and  Thome,  1978]  and  the  other  on  single  particle  motion  [Kivelson  and 
Southwood,  1975;  Ejiri,  1978;  Southwood  and  Kaye,  1979;  among  others].  Both  theories 
assume  convection  due  to  the  cross  tail  electric  field  will  renew  the  plasma  sheet  of 
particles  lost  to  precipitation. 

The  strong  pitch  angle  diffusion  theory  is  based  on  the  assumption  particle  lifetimes 
approach  the  strong  pitch  angle  diffusion  limit  (detailed  in  Chapter  3).  This  facilitates  an 
isotropic  pitch  angle  distribution  where  the  particle  lifetimes  T  approach  a  minimum 
[Kennel  and  Ashour-Abdalla,  1982]. 
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(2-1) 


2Xi/4B 


where  Xi/4B  is  the  quarter  bounce  period  (i.c.  between  the  equator  and  the  atmosphere)  and 
ttL  is  the  width  of  the  equatorial  loss  cone. 

As  particles  convect  toward  the  earth  from  the  magnetotail  few  particles  are  lost 
from  the  flux  tubes.  This  is  a  consequence  of  the  electron  minimum  lifetime  being  much 
longer  than  characteristic  flow  time  Xf  at  great  distances.  In  a  dipole  field, 

Xm  ~  (2-2) 


where  Ec  is  the  electron  energy  and  L  is  the  distance  from  the  earth  in  a  dipole  magnetic 
field  measured  in  Re.  L  is  known  as  the  Mcllwain  parameter  [Mcllwain,  1961].  The 
characteristic  flow  time  is 


where  <}»  is  the  electric  potential  across  the  magnetotail. 

Closer  to  the  earth  the  magnetic  flux  tubes  decrease  in  volume,  thus  increasing  the 
flux  density.  In  (2-2),  Ec  changes  with  L  adiabatically  so  that  finally,  for  characteristic 
electron  energies  at  small  L  (i.c.  ^  7-lORe),  Xm  is  smaller  than  the  characteristic  Xf.  At  this 
point  electrons  rapidly  precipitate  [Kennel  and  Ashour-Abdalla,  1982].  Due  to  their  mass, 
protons  would  be  expected  to  have  longer  lifetimes  and  the  inner  edge  of  the  plasma  sheet 
should  be  closer  to  the  earth.  Yet,  the  electron  inner  edge  is  observed  closer  to  the  earth 
(Newell  and  Meng,  1987].  This  is  explained  by  Ashour-Abdalla  and  Thome  [1978]  who 
suggest  electron  precipitation  in  regions  of  field  aligned  currents  can  create  an  environment 
in  which  ion  cyclotron  waves  can  pitch  angle  scatter  ions.  Though  such  ion  cyclotron 
waves  have  been  documented  by  Gurnett  and  Frank  [1977]  and,  during  substorms,  by 
Cornwall  and  Schulz  [1971],  the  connection  to  electron  precipitation  has  not  been 
confirmed. 

The  second  theory  on  the  formation  of  the  inner  edge  of  the  plasma  sheet  is  based 
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on  the  motion  of  single  particles  in  a  dipole  magnetic  field  and  a  Stern- Volland  cross  tail 
electric  field .  No  initial  diffusion  process  is  assumed.  A  model  of  the  particle  motion  in 
the  equatorial  plane  of  the  magnetosphere  is  given  by  the  electrostatic  potential  [Roederer, 
1970].  As  electrons  convect  toward  the  earth,  zero  energy  particles  follow  electric 
equipotentials  that  exclude  them  from  the  plasmapause.  Higher  energy  electrons  are 
affected  by  eastward  magnetic  drifts  proportional  to  their  energies.  The  larger  the  energy 
the  farther  from  the  earth  the  particle  will  be  deflected.  The  lowest  energy  electrons  will 
convect  closest  to  the  earth  to  form  the  inner  edge  of  the  plasma  sheet.  Since  ions  drift  to 
the  west  they  partially  standoff  the  eastward  corotational  force.  Ions  should  convect  closer 
to  the  earth  than  equivalently  energized  electrons.  This  is  usually  not  observed;  the 
exception  is  at  dusk  [Newell  and  Meng,  1987].  In  this  theory,  precipitation  of  electrons 
into  the  diffuse  aurora  would  require  a  scattering  mechanism  such  as  waves. 

Neither  model  of  magnetotail  convection  is  fully  supported  by  observations.  The 
convection  of  particles  in  the  strong  diffusion,  as  assumed  in  the  first  theory  above,  is  not 
proven,  Fairfield  and  Vifias  [1983]  suggest  that  the  inner  edge  of  the  plasma  sheet  is  a 
function  of  both  theories. 

2.3  RESEARCH  CHRONOLOGY 

A  knowledge  of  the  mean  rate  of  diffusion  is  requisite  in  the  determination  of  the 
primacy  of  ECH  waves  sustaining  the  diffuse  aurora.  In  this  work,  I  investigate  the 
assumption  of  pitch  angle  scattering  at  the  strong  diffusion  limit.  This  is  done  by 
comparing  observed  pitch  angle  distributions  to  those  theoretically  derived  from  this  limit. 
The  ECH  wave  amplitudes  necessary  to  create  the  observed  distributions  in  this  study  arc 
computed.  These  amplitudes  are  then  compared  to  the  ECH  wave  observations  of  Roeder 
and  Koons  [1989]  and  Belmont  et  al.  [1983].  This  thesis  is  another  step  in  a  series  of 
studies  dating  back  to  Kennel  et  al.  [1970]. 
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The  Ogo  5  satellite,  launched  in  March  1968,  observed  several  new  types  of 
magnetospheric  electric  field  emissions  with  wave  frequencies  Ci)  above  the  local  electron 
cyclotron  frequency  fi..  The  most  common  of  these  waves  have  become  known  as 
Electron  Cyclotron  Harmonic  (ECH)  waves  because  of  its  tendency  to  occur  at  odd  half 
multiples  of  ft..  Kennel  et  al.  [1970]  named  ECH  waves  as  a  likely  source  of  pitch  angle 
diffusion  and  turbulent  energization  of  auroral  zone  electrons.  They  cited  the  fact  that  ECH 
waves  with  to  3ft 72  are  observed  near  the  geomagnetic  equates  on  field  lines  anchored  in 
the  auroral  region  and  typically  have  intensities  of  l-lOmV/m.  Lyons  [1974],  using  quasi- 
linear  diffusion  theory,  showed  that  these  waves  should  generally  cause  strong  pitch  angle 
diffusion  of  electrons  of  a  few  tenths  keV  to  a  few  keV.  Lyons  noted  that  the  most  intense 
waves  could  put  electrons  of  energies  up  to  -lOOkeV  on  strong  diffusion.  The  work  of 
Kennel  et  al.  [1970]  and  Lyons  [1974]  is  the  theoretical  basis  to  which  all  observational 
surveys  of  the  plasma  sheet  have  since  been  compared. 

Continued  analysis  of  CXjO  5  data  [Fredricks  and  Scarf,  1973]  and  new  data 
available  from  GEOS  1  and  GEOS  2  (Gough  et  al.,  1979;  Canu,  1982]  afforded  better 
definition  of  the  location,  intensity,  and  frequency  of  occurrence  of  ECH  waves.  It  was 
noted  the  waves  were  more  closely  concentrated  on  the  geomagnetic  equator  than  had 
originally  been  suggested.  The  observations  from  GEOS  indicated  intense  ECH  waves  did 
not  commonly  occur  [Gough  et  al.,  1981]. 

Belmont  et  al.  [1983]  refined  Lyons'  [1974]  theoretical  calculations  by  adapting 
new  information  derived  by  Gough  et  al.  [1979,1981]  and  Canu  [1982].  Wave 
interactions  were  confined  to  within  a  few  degrees  of  the  magnetic  equator,  instead  of 
assuming  interaction  along  the  entire  field  line.  A  minimum  wave  amplitude  required  to 
place  a  resonant  IkeV  electron  on  strong  diffusion  was  determined  to  be  in  excess  of 
2mV/m.  Belmont  et  al.  found  amplitudes  of  2mV/m  were  observed  only  9%  of  the  time 
when  GEOS  was  closest  to  the  magnetic  equator  and  only  1%  of  the  time  at  -3°  off  the 
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equator.  The  paucity  of  ECH  waves  with  adequate  amplitudes  implies  they  are  not  the 
single  (or  even  the  primary)  mechanism  for  putting  IkeV  electrons  on  strong  diffusion 

Fairfield  and  Vinas  [1984]  concluded  pitch  angle  diffusion  does  not  occur  at  the 
strong  diffusion  rate  in  the  plasma  sheet.  They  based  their  conclusions  on  two 
observations.  First,  the  inner  edge  of  the  plasma  sheet  was  too  close  to  the  plamapause. 
Had  electrons  been  on  strong  diffusion  they  would  have  been  lost  sooner  and  the  inner 
edge  of  the  plasma  sheet  would  stand  farther  out  The  second  observation  was  of  "frequent 
and  persistent  occurrences  of  anisotropies  [in  the  electron  distributions]."  Fairfield  and 
Vinas  suggested  that  under  the  influence  of  strong  diffusion  these  anisotropies  would  be 
smoothed.  Thus,  they  support  the  single  particle  theory  for  the  convection  of  electrons. 
This  theory  would  cause  electrons  normal  to  the  magnetic  field  to  be  energized  more  than 
field  aligned  particles.  The  longer  a  particle  distribution  remains  under  the  influence  of  the 
convection  electric  field  in  a  dipole  magnetic  field  the  greater  the  pitch  angle  anisotropy. 

The  electron  distributions  studied  by  Schumaker  et  al.  [1989]  did  not  support 
strong  diffusion  by  ECH  waves  either.  Though  the  distributions  were  isotropic  in  the  loss 
cone,  immediately  outside  the  loss  cone  they  become  anisotropic.  This  type  of  distribution 
would  require  the  diffusion  coefficient  to  be  large  across  the  loss  cone  and  then  decrease 
toward  higher  pitch  angles.  This  is  contrary  to  the  altitude  normalized  diffusion  coefficient 
calculations  of  Lyons  [1974].  Lyons'  calculations  showed  pitch  angle  diffusion 
coefficients  peak  well  away  from  a  =  0°  and  decreases  towards  a  =  0®  (see  Fig.  3-3).  The 
diffusion  coefficient  discrepancy  notwithstanding,  ECH  wave  amplitudes  great  enough  to 
produce  the  distributions  recorded  by  Schumaker  et  al.  have  still  not  been  observed. 

Roeder  and  Koons  [1989]  conducted  a  survey  of  ECH  emissions  using  AMPTE 
and  SCATHA  data  at  ±10®  and  ±5®  of  the  geomagnetic  equatw,  respectively.  They  found 
that  even  in  the  region  of  maximum  ECH  emissions  (i.e.,  0300-0600MLT,  4-8  Re) 
amplitudes  exceeding  0.012mV/m  were  observed  only  60%  of  the  time.  Amplitudes  of 
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0.035raV/m  were  observed  only  25%  of  the  time.  Roeder  and  Koons  concluded  the 
occurrence  and  amplitude  of  ECH  waves  were  too  small  to  create  the  necessary  pitch  angle 
diffusion  to  account  for  the  continuous  diffuse  auroral.  This  conclusion  supported  the 
earlier  work  of  Belmont  et  al. 


10 


PITCH  ANGLE  DIFFUSION 


3.1  ATMOSPHERIC  LOSS  CONE 

Particles  in  a  magnetic  mirror  can  be  characterized  as  precipitating  or  trapped 
depending  on  their  location  in  velocity  space.  A  trapped  particle  mirrors  on  a  field  line 
described  by 

sin  2  00  =  1^  (3-1) 

where  oCq  is  the  equatorial  pitch  angle.  Bo  is  the  magnetic  field  strength  at  the  equator  (the 
weakest  strength  on  a  field  line),  and  is  the  mirror  point  field  strength. 

Precipitating  particles  do  not  mirror  on  their  respective  field  lines  before  entering  the 
earth's  atmosphere  and  collide  with  ionospheric  particles.  The  range  of  pitch  angles  at  the 
equator  who's  mirror  point  is  within  the  earth's  atmosphere  is  defined  as  the  atmospheric 
loss  cone  ai,.  The  height  of  the  earth's  atmosphere  dense  enough  to  significantly  increase 

the  probability  of  collision  is  nominally  100km.  The  atmospheric  loss  cone  is  given  by 

Ol  =  arcsin  (^J^)‘^  (3-2) 

where  Bioo  is  the  magnetic  field  strength  at  100km.  The  loss  cone  width  is  independent  of 
particle  energy,  but  there  is  still  a  particle  species  dependance.  The  larger  cross-section  of 
protons  increases  their  probability  of  collision  at  altitudes  higher  than  electrons  thus  giving 
them  a  wider  loss  cone. 

A  bounce  period  Tb  is  defined  as  the  time  necessary  for  a  particle  to  travel  from  the 
equator,  mirror  in  both  hemispheres,  and  return  to  the  equator.  Particles  with  equatorial 
pitch  angles  Oo  less  than  or  equal  to  the  loss  cone  angle  are  precipitated  into  the  atmosphere 
within  a  quarter  bounce  period  Ti/4b  (i.e.  lost  between  the  equator  and  the  first  mirror 
point).  Particle  lifetimes  Tl  are 


'Cl  =  ''^i/4b  for  do  ^  at 


(3-3) 


Xl  =  oo  for  cXo>aL  (3-4) 

Equation  (3-3)  defines  loss  cone  particles  and  (3-4)  defines  trapped  particles.  Trapped 
particles  in  the  plasma  sheet  are  more  correctly  called  quasi- trapped.  The  are  trapped  until 
such  time  as  they  are  lost  from  the  magnetosphere  or  scattered  into  the  loss  cone. 


3.2  ADIABATIC  INVARIANTS 


Three  periodic  motions  arise  from  trapped  particles;  gyration,  bounce  and  drift. 
The  first  and  fastest  motion  is  a  particle's  gyration  about  a  magnetic  field  line.  The  second 
motion  is  the  periodicity  associated  with  a  particle's  bounce  between  mirror  points.  The 
third  and  slowest  motion  is  the  azimuthal  drift  around  the  earth.  If  the  particle  moves 
slowly  through  variations  in  the  magnetic  field  these  periods  are  conserved. 

Conservation  of  the  magnetic  moment  by  a  particle's  gyromotion  is  called  the  first 


adiabatic  invariant  and  is  defined 

mVj.2  pj2  e  sin^a 

2B  "  2Bm  =  B 


(3-5) 


where  m  is  the  mass,  Vj.  is  the  particle  velocity,  px  is  the  momentum,  both  Vx  and  px  are 
normal  to  the  magnetic  field,  E  is  the  particle  energy,  and  a  is  the  pitch  angle.  The  second 


adiabatic  invariant  is  the  relationship  of  a  particle's  bounce  motion  on  field  lines 

j=  /pii  di 


(3-6) 


where  p  1 1  is  the  parallel  component  of  the  momentum  and  1  is  the  field  line  length  from  the 

equator.  Conservation  of  the  magnetic  flux  in  a  particle's  drift  motion  is  defined  as  the 
third  adiabatic  invariant 
<D  =  JB  dA 
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(3-7) 


where  A  is  the  area  within  the  drift  orbit. 


3.3  PITCH  ANGLE  DIFFUSION 

Processes  that  violate  either  the  first  or  second  adiabatic  invariant  can  scatter,  or 
diffuse,  quasi-trapped  particles  to  lower  equatorial  pitch  angles  and  eventually  into  the  loss 
cone.  Variations  in  electric  or  magnetic  fields  or  interaction  with  magnetospheric  waves 
can  cause  such  diffusion.  Here  we  will  concentrate  on  wave-particle  interactions. 

Diffusion  in  velocity  space  due  to  resonant  waves  interacting  with  particles  is  a 


combination  of  pitch  angle  and  energy  diffusion. 
(Aa)2 


Dan  - 

Dee  = 


2Xb 

(AE)2 


(3-8) 

(3-9) 


E2  2Tb 

where  D  is  the  diffusion  coefficient,  Aa  and  AE  are  the  small  changes  in  pitch  angle  and 
energy,  respectively.  A  simplifying  argument  we  will  assume  is  that  a  particle's  random 
walk  in  energy  diffusion  is  much  slower  than  in  pitch  angle  (Daa»DEE)-  The  reasons  for 
making  this  assumption  is  the  invocation  of  3/2Q.  ECH  waves  in  Section  3-5.  Kennel 
[1969]  states  it  is  not  inappropriate  to  disregard  the  energy  diffusion  and  examine  only 
pitch  angle  diffusion  for  the  small  phase  velocity  of  electrostatic  3/2Q-  ECH  waves. 

The  time  it  takes  a  particle  to  diffuse  across  the  loss  cone  in  a  quarter  bounce  period 
represents  the  separation  of  strong  and  weak  diffusion.  In  weak  diffusion,  particles  will 
fail  to  completely  diffuse  across  the  loss  cone  in  a  quarter  bounce  period  and  are 
precipitated. 

,2 

»  1  (3-10) 


Daa  Xi/4B 

Xd»Xi/4B  (3-11) 

where  to  is  diffusion  time.  The  weaker  Dom,  the  weaker  the  flux  in  the  loss  cone.  In  very 
weak  diffusion  the  loss  cone  is  virtually  empty. 
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In  strong  diffusion,  particles  can  be  diffused  into  and  out  of  the  loss  cone  in  a  quarter 


bounce  period. 

CtL^  , 

(3-12) 

Daa  Ti/4B 

XD  «  Ti/4B 

(3-13) 

The  loss  cone  remains  filled  in  strong  diffusion  since  precipitated  particles  arc  immediately 
replaced. 

Figure  3-1  is  an  enlightening  analogy  to  the  levels  of  pitch  angle  diffusion  as 
illustrated  by  Kennel  et  al.  [1969]  (originally  formalized  by  O'Brian  [1964]).  The  pump 
represents  the  source  of  particles,  the  water  level  in  the  bucket  represents  the  quasi-trapped 
flux  intensity,  the  spout  is  analogous  to  the  loss  cone.  In  3- la.,  the  steady  state  case  is 
shown  where  weak  diffusion  balances  injection  and  precipitation.  The  source  in  this  case 
determines  the  loss  rate  to  atmosphere.  Strong  particle  injections  are  represented  in  3- lb. 
Here  the  pump  overwhelms  the  spout's  removal  capacity.  The  spout  (loss  cone)  regulates 
the  rate  of  loss. 

3.4  PITCH  ANGLE  DISTRIBUTIONS 

The  distribution  of  particle  pitch  angles  is  described  by  the  Fokker-Plank  equation 

if  -  — ^  (Doa  sin  a  :^)  =  S(a,E)  -  L(a3)  (3-14) 

sin  a  9a 

where  f(x,v,t)  is  the  one  particle  (Vlasov)  distribution  function  in  velocity  space,  S(a,E) 
represents  sources,  and  L(aJE)  arc  sinks  in  pitch  angle  and  energy. 

Emulating  Kennel  and  Petschek  [1966],  we  will  solve  for  the  pitch  angle 
distribution  assuming  a  steady  state  diffusion  solution.  Two  preliminary  solutions  will  be 
required;  one  inside  the  loss  cone  and  one  outside.  These  solutions  take  into  account 
differences  in  sources  and  sinks. 


14 


WEAK  DIFFUSION  STRONG  DIFFUSION 

«  KEV  ELECTRONS  AURORAL  ELECTRONS  15  KEV) 


(a)  (b) 


Figure  3-1.  Leaky  bucket,  (a)  In  weak  diffusion,  particles  precipitate  at  the  rate 
they  are  added.  A  steady  state  is  obtained,  (b)  In  strong  diffusion,  the 
precipitation  mechanism  saturates  due  to  the  finite  size  of  the  loss  cone.  [From 
Kennel,  1969] 
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First,  solve  for  the  distribution  in  the  loss  cone  where  we  assume  no  sources 


(S(a,E)  =  0).  The  particle  sink  :n  the  loss  cone  is 
f  f 

L  =  -  =  -^  (3-15) 

'^L  'Cl/4B 

From  (3-14)  steady  state  diffusion  in  the  loss  cone  becomes 

—  ^fDaaa  |^y;^  =  0  (3-16) 

a  daj  Ti/4b 

where  the  law  of  small  angles  is  used  for  a 

Applying  the  assumption  of  Doa  »  Dee.  the  pitch  angle  distribution  in  the  loss  cone 


(3-16) 


is  considered  independent  of  time  and  can  be  expressed 


g(a)~X(E)Io 


(Daa  'Ci/4b)*^ 


(3-17) 


where  lo  is  a  modified  Bessel  function  and  X(E)  is  an  arbitrary  function. 


Outside  the  loss  cone  (a  >aL)  there  are  no  sinks  and  the  steady  state  diffusion 


equation  (3-14)  is 


—  —  fD  sin  a  S(a,E) 
n  a  da  da) 


(3-18) 


where  S(a,E)  is  the  source  flux  of  the  distribution.  In  steady  state,  the  source  flux  will 
equal  the  precipitation  flux.  Kennel  and  Petschek  [1966]  make  the  following  simplifying 
assumption;  particles  are  injected  only  at  flat  pitch  angles  (i.e.  n/2)  and  the  diffusion  flux  is 
constant  at  any  angle  a  <  n/2  and  is  equal  to  S(a,E).  The  source  flux  can  be  broken  into 
components  of  energy  and  pitch  angle,  S(aJS)  =  S*(E)  S(a  -  n/2).  The  flat  pitch  angle 
assumption  requires  S(a)  =  0  everywhere  except  at  a  =  n/2. 


The  equation  of  the  distribution  function  exterior  to  the  loss  cone  which  satisfies 


(3-14)  is 


D*  tan  a  — 


a^  =  ^sin  a  S(a  -  ^)  S*(E)  da 


(3-19) 


=  S*(E) 
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where  D  =  D*/cos  a.  Equation  (3-19)  quantifies  the  assumption  that  the  source  strength  on 
the  right  side  of  the  equation  is  balanced  by  the  diffusion  flux  on  the  left. 


Integrating  to  solve  for  f  yields, 

f(a)  =  In  (sin  a)  +  Zo 

where  Zo  is  a  constant  of  integration. 


(3-20) 


If  f(a)  and  the  diffusion  flux  are  considered  continuous,  then  matching  the  loss 
cone  distribution  (3-17)  with  the  external  distribution  (3-20)  will  eliminate  X(E)  and  Zq. 


The  complete  exterior  distribution,  symmetric  about  a  =  nP.  is  then 


f(a«  =  5^rh(a.).l„f^ 

^  Sin  ttL 


(3-21) 


where  h(ao)  is  the  loss  cone  solution 


h(a)  = 


(D*  Ti/4b)^^  °  l(D*  Ti/4b)>^^ 


(D*  Xi/4b)1/2 


(3-22) 


The  lifetime  for  trapped  particles  is  equal  to  the  total  number  of  trapped  particles 
divided  by  the  panicle  loss  rate.  Recall  from  (3-19),  S*(E)  is  equivalent  to  the  precipitation 
flux  in  the  loss  cone.  Dividing  the  exterior  distribution  (3-21)  by  S*(E)  yields 

%a 


1  f  •  ( X  .  ("sin  a'l't 

'^L  =  n*  Jsin  a  h(ao)  +  In  - 

^  ttL 

1  ("u/  X  .  9  tan  tto 

“  ^h((Xo)  -  In  2  J 


(3-23) 


where  q  is  the  charge. 

Observed  particle  pitch  angle  distributions  in  weak  diffusion  are  highly  anisotropic 
across  the  loss  cone.  Quantitatively,  this  is  bore  out  in  the  asymptotic  representation  of  the 
Bessel  function  where  the  loss  cone  distribution  increases  exponentially  from  the  center  (a 


=  0)  to  the  edge  (a  =(Xo). 
In  weak  diffusion, 
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A  minimum  lifetime  sets  an  upper  limit  on  the  quasi-trapped  particles'  rate  of  removal. 
Note  Tm  is  controlled  by  the  size  of  the  loss  cone  and  is  independent  of  the  strength  of  the 
diffusion  coefficient  D  and  charge  q  (i.e.  particles  can't  be  removed  faster  than 
regardless  of  the  strength  of  the  diffusion). 

The  Kennel  and  Petschek  treatise  on  steady  state  diffusion  is  restrictive  in  that  pitch 
angle  scattering  is  limited  to  only  one  location  on  a  field  line,  and  their  definition  of  Daa  ~ 
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D*/cos  a  is  not  representative  in  the  magnetosphere.  Lyons  et  al.  [1972]  calculated  a  time 
average  D  over  an  electron  bounce  period.  This  was  defined  by  relating  the  local  D(a)  to 
the  equivalent  equatorial  D(ao).  The  bounce  averaged  equatorial  pitch  angle  diffusion 


coefficient  is 


Tb 


(3-32) 


Equation  (3-32)  can  be  transformed  over  magnetic  latitude  X  to  produce 
.  X  1  r  .  cos  a 


D(E,L,ao)  = 


s(ao) 


D(a)  cos^XdX 

COS^Oo 


(3-33) 


where  Xm  is  the  mirror  latitude,  and  s(ao)  is  a  correction  for  the  variation  in  Tb  for  a  given 
a  (s(<Xo)=  1.38  -  .032(sin  Oo  +  sin*/2(x^j)  [Schulz  and  Lanzerotti,  1973]. 

Next  f(x,v,t)  is  converted  from  velocity  space  to  phase  space  f(E,L,ao,t)  by 
averaging  (3-14)  over  an  electron  bounce  in  a  dipole  field. 

^  ^ - ^  s(ao)  sin  (oto)  cos  D(ao)  (3-34) 

s(ao)  sin  ao  cos  Cq  dOo  dOo  ^atm 

where  Tatm  =  Xi/2B  (half  bounce  period)  for  a  <  Oq  and  Tami  =  «» for  a  >  at- 


The  steady  state  model  is  rarely  seen  in  the  magnetosphere.  More  common  are 
electron  injections  on  the  order  of  minutes  followed  by  several  hours  of  decay.  Roberts 
[1969]  studied  artificial  injections  by  a  nuclear  explosions  and  natural  injections  and  found 
they  both  soon  decayed  to  a  normal  mode.  Lyons  used  this  fact  to  break  the  equatorial 
pitch  angle  distribution  into 

fo(t.ao)  =  F(t)g((Xo)  (3-35) 

where  F(t)  denotes  the  equatorial  density  of  particles  for  a  given  energy  and  gfOo)  denotes 


the  shape  of  the  distribution. 
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Substituting  (3-35)  into  (3-34)  produces 
f  —  +  t  ^  1  g(«o)  s(ao)  sin  cto  cos  Oo 

V  Mtm  y 

=  s(ao)  sin  Oo  cos  Oo  D(ao)  ^ 
doc^  dcio 


(3-36) 


Lyons  converted  (3-36)  to  integral  form  by  carrying  out  the  integration  twice 

between  0  and  n/2,  specifying  =  0  at  Oo  =  0®  and  Oo  =  90°.  The  resulting  equatorial 

3ao 

pitch  angle  distribution  g(ao)  takes  the  form 

tto 

(* _ I _  . 

g(0)  To  J  g(a')  s(a')  sin  a’  cos  a’  D(a’) 


/  a”  a' 

Jg(a)  sin  a  cos  a  da  -  q  Jg(a)s(a)  sin  a  cos  a  da') 


(3-37) 


where 


a"  =  a'  for  a'  <  at 
a"  =  aL  for  a'  >  at 


Jg(ao)  sin  a 


cos  ao  da. 


(3-38) 


Jg(a« 


)  s(ao)  sin  ao  cos  ao  da^ 


Thus,  the  equatorial  pitch  angle  distribution  can  be  determined  for  given  values  of  DfOo). 


The  precipitation  lifednoe. 


1  aFvi 


/r  /"I  o  n 

^p=-(F5r) 


is  related  to  gfOp)  by 
T 

Vn  — 


(3-39) 


(3-40) 
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L-8 


'ii{.o)a 
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Equotonol  Pitch -Angle 


where  To  =  4LRe/v  (v  is  particle  velcx:ity)  and  T|  is  the  probability  of  a  particle  being  in  the 
loss  cone.  T]  is  equal  to  the  integral  of  the  loss  cone  over  the  integral  of  the  entire 


Add'mg  in  the  field  line  correction  2s(ao)  for  a  dipole  field 

q  =  2s(ao)  sin^at  (3-42) 

(3-43) 

n 


Lyons  [1973]  went  on  to  calculate  numerical  solutions  for  (3-37)  to  produce 

theoretical  equatorial  pitch  angle  distributions.  He  assumed  no  sources  and  a  constant  pitch 

angle  distribution  shape  as  the  flux  decayed.  Values  of  D(0“)  and  D(90°)  were  chosen  with 

D(oto)  allowed  to  vary  linearly  between  tto  =  0®  and  a®  =  90®.  Note  that  the  delineation 

OtL^ 

between  weak  diffusion  (3-10)  and  strong  diffusion  (3-12)  is  approximately  — - =1; 

Daa  Ti/4B 

cil^ 

where  Daa  = - “  By  normalizing  D(0®)  to  Xn,  *,  solutions  to  (3-37)  can  be 

Ti/4B 

separated  into  diffusion  regimes  of  very  weak  (D(0®)Tm  =  0.01),  weak  (0.1),  moderate  (1), 
and  strong  (10).  Plotted  in  Figure  3-2  are  Lyon’s  solutions  to  (3-37)  for  L  =  2,  4,  and  8. 


3.5  ELECTROSTATIC  CYCLOTRON  HARMONIC  WAVES  AND  DIFFUSION 

The  diffusion  process  presented  thus  far  has  been  general  and  no  scattering 
mechanism  has  been  specified.  Now  we  apply  the  Kennel  et  al.  [1970]  suggestion  that 
ECH  waves  are  responsible  for  pitch  angle  scattering  in  the  near  earth  plasma  sheet.  Lyons 
[1974]  quantified  the  Kennel  et  al.  theory  for  3/20.  ECH  waves  using  quasi-linear 
diffusion.  Kennel  and  Englemann  [1966]  presented  a  general  quasi-linear  diffusion 
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equation  acted  upon  by  a  specific  wave  distribution.  That  equation  for  a  space  plasma  in 


the  presence  of  ECH  waves  is 


af 

g^  =  hmv. 


±nft  a 


(3-2703  Vlk|2 

Jn^(kiVi/0)  / 
vk  -  k||V|i  ±  (-nft)  { 


lk|2  ^  vjL  avi  ~  < 

/  ±nQ  a  1,  3  \ 

vjL  avj.  ~  ’*  5^  j 


a  N 
5^  J* 


(3-44) 


where  f  is  the  spatially  uniform  zeroth-oider  distribution  function  in  phase  space,  m  is  the 
plasma  species  rest  mass,  q  is  the  charge,  and  ±  refers  to  the  sign  of  the  charge.  V  is  the 
plasma  volume,  vk  =  cok  -t-  yii  is  the  complex  wave  frequencies  as  a  function  of  the  wave 
number  k,  Ek  is  the  electric  field  at  each  k,  Q  =  (qB/mc)  is  the  gyrofrequency,  c  is  the 
speed  of  light,  and  Jn  is  a  Bessel  function  of  order  n. 

Equation  (3-44)  can  be  transformed  from  f(vj.,V||)  to  f(v,a) 

17  .  /Til  ^  1  ~  /  t\  •  ^  .  r\  » 


=  V.(D'.Vf)=— sinaf  D„a*  |^  +  D„v’ 

V  sin  a  aa  ''  aa  °  J 

la  2/'r» 

5^J 


(3-45) 


where  Dw  is  the  diffusion  coefficient  in  terms  of  speed.  We  are  specifically  concerned 


with  the  diffusion  coefficient  for  electron  interactions  with  electrostatic  3/2  ft.  waves. 


Doa  —  ^  Doatt 


(3-46) 


The  pitch  angle  diffusion  coefficient  is  constrained  by-*^'^^^  ^  ^  and 


Doa"  =  0  eveiywhere  else. 


Averaging  Daa  over  a  particle  bounce  period,  as  discussed  previously  in  (3-32), 
Lyons  [1974]  was  able  to  calculate  numerical  solutions  for  a  normalized  pitch  angle 


coefficient. 


^aa  (Cto»K/Kth)  —  ir?  *i** 


lEwlo^ 


(3-47) 


where  K  is  the  particle  energy,  Kth  is  the  thermal  energy,  lEwl©  is  the  equatorial  wave 
amplitude,  and  <Daa>  is  the  bounce  averaged  diffusion  coefficient.  The  results  of 
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normalizing  the  diffusion  coefficient  to  Kuj/lEwlo^  are  shown  in  Figure  3-3.  He  assumed 
Kthi=  Kihii  and  that  K  remain  constant  The  minimum  3/20.  wave  amplitudes  required  to 
produce  strong  diffusion  for  various  values  of  K/Kth  at  L=7  are  plotted  in  Figure  3-4. 
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^aa  L-7 


Figure  3-3.  Normalized  coefficients  for  diffusion  in  pitch  angle  (Daa)  are  plotted 
as  a  function  of  equatorial  pitch  angle  at  L  =  7  for  representative  values  of  the 
electron  energy  normalized  to  the  thermal  energy  of  the  warm  electrons.  [From 
Lyons,  1974] 
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Figvirc  3*4.  The  minimum  equatorial  wave  amplitude  required  for  strong  pitch 
angle  diffusion  at  L  =  7  is  shown  as  a  function  of  electron  energy  for  dermal 
energies  of  0.1, 1,  and  lOkeV.  [From  Lyons,  1974) 
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CHAPTER  4. 

INSTRUMENTS 


4.1  SCATHA  (P78-2)  SATELLITE 

The  Spacecraft  Charging  AT  High  Altitude  (SCATHA)  satellite,  mission  P78-2,  is 
a  joint  United  States  Air  Force/NASA  vehicle.  The  primary  objective  of  the  satellite  is  to 
provide  environmental  and  engineering  information  on  the  charging  of  spacecraft  surfaces. 
A  secondary  objective  is  the  collection  of  scientific  data  and  is  not  necessarily  related  to 
spacecraft  charging.  A  payload  of  nine  experiments  were  provided  by  various  members  of 
the  scientific  community  [Stevens  and  Vampola,  1978]. 

SCATHA  was  launched  aboard  a  Delta  rocket  on  30  January  1979.  The  satellite 
reached  final  orbit  on  2  Febmary.  P78*2  is  in  a  near  synchronous  orbit  generally  within  8 
degrees  of  the  geographic  equator.  The  orbit  is  elliptical  with  apogee  at  ~  57000km  (R  ~ 
8.9  Re)  and  perigee  at  ~  33200kin  (R  ~  5.2  Re)  geocentric.  The  orbit  is  inclined  ~  7.8 
degrees  from  the  earth's  equatorial  plane.  SCATHA  has  an  orbital  period  of  ~  23.6hrs. 
This  period  results  in  an  eastward  longitudal  drift  of  ~  5.3  degrees  per  day  (see  Fig.  4-1). 

The  P78-2  satellite  body  is  cylindrical  in  shape  with  both  a  length  and  diameter  of 
~  1.75m.  The  body  is  divided  by  a  central  (bellyband)  portion  which  contains  most  of  the 
experimental  payload.  Instruments  sensitive  to  spacecraft  influences  are  isolated  on  booms 
projecting  firom  the  bellyband.  Shown  mounted  fore  and  aft  of  the  bellyband  in  Figure  4-2 
are  solar  cells.  SCATHA  design  and  construction  is  further  detailed  in  Stevens  and 
Vampola  [1978]  and  FenneU  [1982]. 

SCATHA  spins  about  its  axis  at  a  rate  of  ~  Irpm.  During  the  period  of  interest 
(1979-82),  the  spin  axis  was  in  the  plane  of  orbit  directed  towards  dawn.  (In  1986,  the 
axis  was  lifted  to  its  present  position  perpendicular  to  the  orbital  plane,  in  an  effort  to 
extend  the  life  of  the  satellite.)  The  spin  axis  had  to  be  maintained  approximately  normal  to 
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MOTION  OF  P78-2  RELATIVE  TO 
A  SYNCHRONOUS  SATELLITE 


SYNCHRONOUS 


DRIFT  ~5.3°l DAY 


Figure  4-1.  SCATHA  orbital  motion.  The  SCATHA  satellite  has  a  23.6  hour 
orbital  period  resulting  in  a  drift  of  ~5.3®  per  day.  The  subsatellite  ground  track  is 
~40®  wide.  [From  Fennell,  1982] 
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CATHA  satellite  and  payload 
- 1  view.  [From  Stevens  and  Vai 


the  satellite-sun  line  to  maximize  the  efficiency  of  the  solar  cells.  This  orbital  geometry 
required  near  weekly  corrections  to  the  satellite’s  orientation  as  the  earth  revolved  about  the 
sun.  This  geometry  generally  has  the  spin  axis  perpendicular  to  the  magnetic  field  direction 
near  local  midnight 

4.2  SC2-3  ELECTROSTATIC  ANALYZER 

The  data  for  this  paper  was  collected  by  the  SC2-3  Electrostatic  Analyzer  (ESA) 
aboard  the  SCATHA  satellite.  The  SC2-3  ESA  is  one  of  a  troika  of  Electrostatic  Analyzers 
(i.e.  SC2-1,  SC2-2,  and  SC2-3).  The  primary  objective  of  the  ESA  is  to  provide  electron 
and  ion  distribution  functions,  between  a  few  eV  to  ~20keV,  at  three  positions  in  the 
spacecraft  plasma  sheath.  The  SC2-3  is  located  on  the  bellyband  with  a  360°  field  of  view 
perpendicular  to  the  spin  axis. 

The  ESA  consists  of  concentric  cylindrical  plates  oppositely  biased  to  deflect 
particles  of  the  appropriate  energy  per  charge.  The  plates  are  stepped  through  a  cycle  of 
biases.  The  center  plate  in  Figure  4-3  is  biased  from  ~  0  to  -t-16(X)V,  while  the  two  outside 
plates  are  simultaneously  biased  from  ~  0  to  -1600V.  The  plates  are  serrated  to  eliminate 
particles  that  impact  the  sides  and  to  reduce  scattered  light  response.  Particles  succeeding 
in  traveling  the  127°  curved  path  are  detected  in  a  Spiraltron  channel  election  multiplier. 

Since  electron  data  was  used  exclusively  in  the  thesis,  no  further  discussion  of  the 
ESA  ion  measurement  capability  is  presented  here.  The  electron  geometric  factor  is  ~1.7  x 
10"^  cm^  ster.  The  electron  energy  resolution  is  AE/E  ~  0.09. 

The  electron  spectrum  is  divided  into  21  channels.  The  number  of  channels 
operating  is  programmable.  Each  channel  collects  data  for  101 .5msec  after  a  ~23msec  bias 
stabilization  period.  A  complete  21  channel  spectrum  could  be  scanned  in  3  seconds.  The 
channels  are  subdivided  into  three  programs  as  shown  in  Table  4-1.  One  program  can  be 
scanned  in  Isec.  The  usual  procedure  during  the  period  of  interest  was  to  only  scan  two 
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vibj 9™  Analyzer,  SC2.  (From  Slovens 


SC2-3  ELECTROSTATIC  ANALYZER 


STEP 

ELECTRON  ENERGY 
(eV/charge) 

PROGRAM  1 

0 

187 

1 

446 

2 

1090 

3 

2580 

• 

4 

4520 

5 

10950 

6 

19400 

7 

RETURN TO ZERO 

PROGRAM  2 

0 

BKG 

1 

87 

2 

316 

3 

815 

4 

1940 

5 

5900 

6 

14400 

7 

RETURN TO ZERO 

PROGRAM  3 

0 

BKG 

1 

17 

2 

40 

3 

612 

4 

1440 

5 

3410 

6 

8200 

7 

RETURN  TO  ZERO 

Tabic  4-1 .  SC2  Electrostatic  Analyzer  Electron  Energies 


I 


programs  at  a  time.  This  meant  14  channels  were  sampled  30  times  per  minute. 

The  ESA  full  width  10%  maximum  electron  angular  response  was  ~9®  parallel  and 
~7°  perpendicular  to  the  spin  axis.  The  full  width  half  power  points  were  -5.5°  parallel  and 
-5.8°  perpendicular.  The  instrument  response  is  represented  as  a  product  of  a  step 
function  in  the  parallel  and  a  Gaussian  function  in  perpendicular  directions,  shown  fitted  to 
the  data  in  Rgure  4-4. 

The  response  function  becomes  particularly  significant  when  trying  to  observe 
electron  pitch  angle  distribution  features  smaller  than  the  angular  response  (e.g.  the  loss 
cone).  Consider  the  geometry  of  Figure  4-5,  where  y  is  the  angular  response,  ocu  is  the 
detector  pointing  angle,  B  is  the  magnetic  field,  and  4  is  the  phase  angle.  For  pitch  angles 
smaller  than  the  angular  response  when  the  detector  is  field  aligned  (Fig.  4-5a),  ^  is 
completely  contained  within  the  angular  response  y.  This  is  contrasted  with  the  pitch 
angles  where  the  detector  is  not  field  aligned  (Fig.  4-5b).  The  result  of  integration  over  ^ 
on  small  pitch  angles  is  a  large  variation  in  the  effective  instrument  response  function  as 
seen  in  Figure  4-6.  At  look  angles  Oj  greater  than  y,  the  effective  instrument  response  is 
constant  This  instrument  response  function  must  be  removed  to  gain  a  true  measure  of  the 
pitch  angle  distribution.  An  algorithm  to  remove  the  response  function  is  discussed  in  the 
Chapter  5. 

Determination  of  field  alignment  is  ascertained  by  a  triaxial  fluxgate  magnetometer 
(SCI  1).  The  magnetometer  is  located  at  the  end  of  a  4m  boom  (see  Fig.  4-2).  The  axes  of 
the  sensor  are  accurately  aligned  with  the  vehicle.  The  SCI  1  samples  four  times  a  second 
with  a  resolution  of  O.Sy.  The  data  was  fit  to  the  Olsen-Pfitzer  [1974]  magnetic  field 
model. 
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PARALLEL 


PERPENDICULAR 


Figure  4-4.  Representation  of  Electrostatic  Analyzer  as  a  product  of  a  step  function 
pai^lel  to  the  spin  plane  of  the  spacecraft  and  a  Gaussian  function  perpendicular  to 
the  spin  plane.  [From  Roeder  and  Gomey,  1986] 


Z  B 


Figure  4-5.  Geometry  of  the  ESA  pointing  angle  Od  and  its  effect  on  the  pitch  angle 
response  function,  (a)  Pitch  angle  a  is  smaller  than  the  angular  response  y  and 
will  be  completely  contained  within  y  when  the  pointing  angle  is  field  aligned  (i.e. 

0°).  Integration  of  phase  angles  ^  at  0®  goes  to  0,  causing  a  significant  variation  in 
instrument  response  (see  Fig.  4-6).  (b)  At  large  pointing  angles,  the  variation  in 
phase  angle  integration  is  smaller.  (After  Roeder  and  Gomey,  1986] 
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Figure  4-6.  Pitch  angle  response  function  of  the  ESA  for  four  detector  pointing 
angles.  Note  the  large  variation  in  response  at  ct^  =  0®.  [From  Roeder  and  Gomey, 
1987J 
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CHAPTER? 

DATA  ANALYSIS  AND  OBSERVATIONS 


5.1  ANISOTROPY 

The  first  criteria  for  selection  of  data  points  was  the  observation  of  the  loss  cone. 
Observation  of  electrons  in  and  near  the  loss  cone  is  essential  for  any  discussion  of  pitch 
angle  diffusion.  To  ensure  the  ESA  response  included  the  loss  cone  (-2.5°  at  geostationary 
orbit)  data  points  were  selected  only  when  the  spin  axis  (S)  was  perpendicular  to  the  local 
magnetic  vector  B  (90  ±  5°).  This  relationship  between  S  and  6  was  required  to  exists  for 
a  minimum  of  15  minutes  to  allow  adequate  time  for  the  collection  of  a  statistically  reliable 
sample.  When  the  relationship  remained  for  longer  periods,  a  15  minute  window  centered 
on  the  mid-point  was  chosen.  A  total  of  617  nightside  points  were  found  between 
February  1979  and  December  1982. 

Seven  electron  channels  were  selected  for  examination;  0.187, 0.446, 1.09,  2.58, 
4.52,  10.95,  and  19.4keV.  This  energy  range  provides  excellent  coverage  of  the  electrons 
primarily  responsible  for  producing  the  diffuse  aurora.  These  channels  are  also  usually 
always  available  in  the  ESA  data  record. 

ESA  plots  of  parallel  "J||"  and  perpendicular  electron  counts  for  the  period 

when  5  was  nearly  perpendicular  to  B  were  reviewed.  The  plots  consist  of  an  average 
plotted  at  each  minute  for  Jj^(70- 1 1 0°),  J||(0-30°),  and  J||(  1 50- 1 80°).  The  mean  includes  the 

conjugate  angles  (181-360°).  Figure  5-1  is  a  sample  of  ESA  data  where  is  plotted  as  a 

dot  and  J||(0-30°,  150-180°)  as  larger  symbols.  Electron  counts  per  second  arc  plotted  on  a 

logarithmic  scale  which  prevented  simple  'eyeball'  averaging  for  a  15  minute  period. 
Instead,  the  value  of  the  sample  counts  are  median  values  for  tightly  grouped  data  points 
and  selectively  chosen  values  that  best  represent  the  period  for  widely  scattered  data.  These 
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values  were  used  to  determine  the  ratio  R  of  perpendicular  to  parallel  electron  counts  in  the 
pitch  angle  distribution. 

R  =  J^/J„  ±  o  (5-1) 

The  uncertainty  a  was  determined  assuming  a  Poisson  distribution.  Thus  the  total  number 
of  counts  occurring  within  the  15  minute  window  must  be  calculated.  This  total  count  J15 
was  computed  for  both  the  parallel  and  perpendicular  direction. 

JlSxMl  “  '^*±•11  ^  P  (5—2) 

where  is  the  electron  count  per  second  selected  to  represent  the  15  minute  data 
sample,  t^  is  the  actual  sampling  time  (i.e.  0. 1  sec),  t^  is  the  time  of  the  observation  (i.e. 

15min),  and  p  is  the  average  number  of  samples  per  minute.  The  number  of  samples  is  a 

function  of  the  number  of  energy  channels  programmed  to  be  collected  and  the  number  of 
degrees  available  for  sampling  in  each  direction  (i.e.  J^=  A40®  versus  J|  =  A30°).  The  most 

common  program,  used  during  72%  of  the  observations,  collected  a  count  for  a  specific 

channel  once  every  two  seconds.  Lesser  used  programs  collected  counts  once  every  three 

seconds  (26%)  and  once  every  one  second  (2%).  In  the  case  of  a  sample  every  two 
seconds  ,  would  be  sampled  an  average  of  6.7  times  per  minute  (i.e.  6.7  =  lsample/2sec 


X  (40°  X  2)/6°sec  i).  iTie  uncertainty  is  taken  for  the  ratio  of  total  number  of  counts 


aJM  l± 


(5-3) 


The  quoted  o  in  this  study  is  the  average  of  the  largest  and  smallest  uncertainty. 


(5-4) 


A  trend/stability  index  was  determined  at  each  data  point  for  every  channel.  The 

index  noted  trends  towards  increasing  or  decreasing  electron  counts  as  well  as  the  stability 
of  the  trend  for  both  and  J^,  individually.  This  index  was  required  for  further  analysis 

when  relatively  stable  periods  would  be  required.  The  pitch  angle  anisotropy  (Jj^^J,,). 
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the  distribution  was  also  indexed  in  a  similar  manner. 

The  decision  as  to  whether  the  observation  was  taken  within  the  plasma  sheet  was 
subjective  and  is  a  source  of  uncertainty.  This  uncertainty  increases  in  quiet  time 
observations  taken  at  the  satellite's  lowest  altitudes  (i.e.  L<6),  particularly  in  the  pre- 
midnight  bins.  The  points  in  these  bins  are  not  always  within  the  mean  plasma  sheet  inner 
edge  extrapolated  by  Gussenhoven  et  al.  [1981].  Rapid  decreases  in  the  electron  counts 
were  considered  indications  of  the  satellite  moving  out  of  the  plasma  sheet  towards  the 
plasmasphere.  Extremely  low  or  erratic  count  rates  were  eliminated. 

Approximately  3400  pitch  angle  ratios  for  7  energy  channels  at  493  different  points 
in  the  plasma  sheet  were  computed  and  considered  valid.  The  spatial  distribution  of  data 
points  on  the  nightside  is  shown  in  Figure  5-2.  The  bulk  of  the  distribution  is  in  between 
2200MLT  and  0200MLT  but  there  is  adequate  and  nearly  uniform  coverage  in  the  pre¬ 
midnight  and  post-midnight  local  times.  The  nightside  coverage  ranged  in  altitude  from 
~5.3L  to  ~8.8L. 

Data  points  were  binned  into  various  ranges  of  local  time,  geocentric  radius  L,  and 
geomagnetic  activity  levels.  Local  time  bins  included  midnight  (2200-0200MLT),  pre¬ 
midnight  (1800-2200MLT),  and  post-midnight  (0200-0600MLT).  L  bins  were  divided 
into  three  ranges;  L  <  6,  6  <  L  <  7,  L  ^  7.  When  combining  the  criteria  for  local  time  and 
L  value  nine  bins  are  formed  as  shown  in  Figure  5-3. 

The  three  hourly  Kp  index  was  selected  for  determination  of  geomagnetic  activity 
binning  because  of  its  availability  and  relation  to  conjunctive  boundaries  of  the  diffuse 
aurora  and  the  inner  edge  of  the  plasma  sheet  [Gussenhoven  et  al.,  1981;  1983]. 
Observations  occurring  during  periods  of  Kp  ^  3  were  classified  as  quiet  and  those  with 
Kp>3  as  active.  Since  the  near  earth  plasma  sheet  may  not  return  to  pre-substorm  electron 
densities  for  several  hours  after  an  injection  [Roberts,  1969;Hoffman  and  Burch,  1973], 
these  definitions  were  expanded.  To  account  for  possibly  elevated  post-injection  electron 
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Figure  5-2.  Spatial  distribution  of  data  points,  (a)  Seven  energy  channels  at  493 
different  points  in  the  plasma  sheet  were  sampled.  Observations  occurring  during 
periods  of  Kp  ^  3  were  classified  as  quiet  (b)  and  those  with  Kp  >  3  as  active  (c). 
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levels,  quiet  time  observations  occurring  within  the  first  two  hours  of  a  3-hour  Kp  period 
required  further  examination.  In  these  cases  if  the  current  Kp  index  was  quiet  but  the 
preceding  Kp  period  was  active,  the  observation  was  binned  as  active.  This  insured  that 
any  observation  classified  as  quiet  was  taken  at  least  two  hours  after  a  substorm  injection. 
This  biased  the  less  frequently  occurring  active  observations  toward  quiet  levels  since  a 
data  point  may  occur  as  much  as  five  hours  after  an  injection  and  still  be  classified  as 
active. 

Table  5-1  contains  statistics  for  the  bins  in  Figure  5-3  for  all  Kp  values  and  Table  5- 

2  for  Kp<3.  The  mean  bin  pitch  angle  ratio  Rbin  is  given  as  well  as  the  average  uncertainty 

Obin-  Mean  bin  ratios  were  calculated  as 
ZR- 

Rbin  =~^  ±  ctbin  (5-5) 

where  n  is  the  number  of  observations  and  the  uncertainty  is  defined 
Zoi 

Obin=—  (5-6) 

Isotropy  is  defined  as  the  percentage  of  observed  bin  samples  where  R  is  within 
one  standard  deviation  of  R=l.  This  standard  deviation  is  determined  from  the  average  o 
of  the  R=1  samples  from  that  bin.  In  cases  where  there  were  no  isotropic  samples  present, 
values  of  equal  energy  from  adjacent  bins  were  used.  Of  the  375  quiet  time  samples 
isotropy  varies  from  0  to  30%  with  an  average  of  8  ±  7.7%.  The  addition  of  1 18  active 
time  samples  increased  the  highest  percentage  of  isotropy  to  33%  and  increased  the  overall 
average  isotropy  to  9  ±  8. 1%. 

The  mean  quiet  and  combined  (quiet  and  active)  pitch  angle  distributions  are  all 
anisotropic  (Jx>J||).  The  average  anisotropies  arc  plotted  in  Figure  5-4  and  show  several 

general  trends.  The  weakest  anisotropies  (more  isotropic)  are  at  the  largest  L  values  (i.e.  L 
^  7).  The  anisotropy  change  is  greatest  between  the  L>7  and  6<  L  <7  regions  at  all  local 
times.  Generally,  the  largest  anisotropies  are  in  the  post-midnight  sector.  These  cases  of 
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AVERAGES:  anisotropy  (Rjin-Ji/Jh),  uncertainly  (obin).  and  %  isotropy 
for  combined  quia;  and  active  times 
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TABLE  $-2  BIN  AVERAGES:  anisotropy  (Rbin-Jx/J||).  uncortainty  (obin).  «nd  %  Isotropy 

for  quiot  tint*  (KPs3) 
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increasing  anisotropy  might  be  expected  when  an  electron's  trajectory,  down  the 
magnetotail  and  eastward  around  the  earth,  is  considered.  Electron  trajectories  computed 
by  Ejiri  [1978]  show  particle  drift  in  the  post-midnight  sector  slows  down  dramatically.  It 
is  hypothesized  by  Fairfield  and  Vifias  [1983]  that  those  particle  distributions  farthest  along 
the  drift  trajectory  have  been  acted  upon  by  anisotropizing  forces  for  a  longer  period  of 
time.  But  this  simple  hypothesis  fails  in  two  other  areas.  First,  the  anisotropy  increases 
as  the  particles  move  from  distant  magnetotail  radially  inward  into  the  middle  bins  (6^  L 
>7)  then,  in  some  cases,  unexpectedly  decreases  as  the  particles  continue  earthward  to 
L<6.  Anisotropy  decreases  at  all  energies  in  all  the  pre-midnight  sector  and  in  half  of  the 
energies  in  the  post-midnight  sector  between  6^  L  <7  and  L<6.  Despite  the  decrease  in 
anisotropy  it  is  greater  in  the  innermost  bins  than  in  the  L>7  bins.  Second,  the  anisotropy 
is  greater  in  the  pre-midnight  bins  than  in  the  azimuthally  downstream  midnight  bins.  This 
latter  contradiction  may  be  evidence  of  plasmaspheric  particles  erroneously  identified  as 
being  in  the  plasma  sheet  and  adding  to  the  anisotropy. 

During  active  times  isotropy  is  more  prevalent.  So  as  expected,  Rbin  at  Kp^3  is 
generally  more  anisotropic  than  Rbin  at  combined  active  and  quiet  Kp.  This  is  particularly 
true  at  the  lower  energies.  The  higher  energies  (10.95  and  19.4keV)  at  quiet  times  show  a 
mix  of  increasing,  decreasing,  and  stable  anisotropies  relative  to  the  same  energies  in  the 
combined  Kp  bins. 

The  pre-midnight  bins  at  6:S  L  <  7  for  the  two  lowest  energies  (187  and  446eV) 
appear  anomalously  more  anisotropic.  This  is  due  to  a  quarter  of  the  pitch  angle  ratios  in 
these  bins  possessing  anisotropies  greater  than  7.  These  values  appear  to  be  valid  but  can 
not  be  immediately  attributed  to  a  specific  phenomena.  A  follow-on  investigation  will  be 
required  to  explain  their  appearance.  When  these  high  values  are  removed,  the  average 
anisotropy  falls  to  levels  more  in  line  with  those  observed  at  higher  energies. 

An  attempt  to  find  a  functional  form  to  represent  the  various  levels  of  anisotropies 
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was  tried.  Fits  to  ySin^a  and  xSin"kx  ±  ySin"a  were  initially  promising  but  limited  time 
to  complete  this  thesis  precluded  a  detailed  study  of  such  representations. 

5.2  DIFFUSION 

To  examine  the  pitch  angle  distributions  in  the  loss  cone  and  determine  the  extent  of 
diffusion,  actual  IS  minute  data  records  had  to  be  examined.  The  data  records  consists  of 
the  time  in  seconds,  average  pitch  angle  for  that  second  and  the  counts  per  0.1  sec  per 
channel.  The  process  of  retrieving  specific  15  minute  periods  of  data  from  the  archived 
tapes  was  time  consuming.  If  for  some  reason  a  data  set  was  unavailable,  incomplete,  or 
otherwise  unacceptable  only  a  few  could  be  replaced  in  a  timely  manner.  Acting  within  this 
constraint  the  19.4keV  energy  channel  was  initially  selected  for  detailed  analysis.  This 
choice  was  an  effort  to  avoid  a  local  maximum  in  the  distribution's  loss  cone  (i.e.  an  anti¬ 
loss  cone)  created  by  secondary  and  backscattered  electrons.  This  phenomenon  is  often 
present  in  energies  up  to  a  few  keV  [Mcllwain,  1975].  Selection  of  representative  sample 
distributions  from  each  quiet  time  bin  at  19.4keV  were  based  on  three  criteria.  First,  the 
particle  counts  must  have  had  a  monotonic  trend  and  stability.  This  eliminated  the  need  to 

adjust  the  data  for  temporal  changes  during  the  15  minute  sampling  period.  Second,  the 
trend/stability  of  the  anisotropy  also  had  to  be  monotonic.  Finally,  ratios  had  to  be 

within  10%  of  the  mean  anisotropy  of  their  respective  bin. 

To  obtain  the  true  electron  distribution,  a  deconvolution  technique  was  applied  to 
the  data  to  remove  the  effects  of  the  detector  response  function.  This  enabled  features 
smaller  than  the  detector's  angular  response  width,  such  as  the  loss  cone,  to  be  discerned. 
The  technique  used  was  specifically  developed  for  the  SCATHA  ESA  data  and  was 
modeled  after  maximum-entropy  schemes  developed  for  image  restoration  and  X-ray 
spectral  deconvolution. 

The  statistical  technique  assumes  the  individual  observations  of  electron  flux  are 
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related  to  the  actual  pitch  angle  distribution  through  a  convolution  function.  The 
convolution  function  describes  the  response  of  the  particle  detector  as  a  function  of  pitch 
angle.  The  deconvolved  distribution  can  be  obtained  by  maximizing  the  entropy  of  the 
electron  distribution  within  the  constraints  of  the  convolution  function  and  a  Chi-squared 
error  function.  The  algorithm  is  fully  described  by  Roeder  and  Gomey  [1986]. 

To  facilitate  the  deconvolution  program,  data  from  the  individual  0.1  sec  count 
periods  between  91-180®  were  'folded  over'  onto  the  observations  between  0-90°.  Then 


the  counts  C  were  averaged  over  2°  bins 
Cl”  =-^±<Jc^ 


(5-7) 


The  uncertainty  ac2»  was  determined  by 


(5-8) 


The  solid  line  in  Figure  5-5  is  an  example  of  a  reconstruction  of  the  pitch  angle 
distribution  in  1°  increments  produced  by  the  deconvolution  algorithm.  The  counts  are 
plotted  as  crossed  error  bars.  The  vertical  error  bar  represents  ac2»  attd  the  horizontal  bar 


is  the  instrument's  full  width  response.  The  slight  undulations  in  the  plot  are  a 
manifestation  of  the  2®  averaging  and  are  not  considered  real.  In  every  case  recomputation 
of  the  pitch  angle  ratio  from  the  reconstructed  distribution  produces  a  less  anisotropic 
distribution  than  obtained  from  using  only  the  Jj^/J||  ratio.  Anisotropies  of  19.4keV 


electrons  decreased  on  average  27  ±  9%. 

The  extent  of  diffusion  was  determined  by  comparison  of  the  deconvolved  pitch 
angle  distribution  to  the  theoretical  distribution  computed  by  Lyons  [1973]  (see  Fig.  3-2). 
Lyons'  predicted  equatorial  pitch  angle  distributions  were  digitized  for  this  purpose.  The 
theoretical  diffusion  coefficients  plots  were  then  normalized  to  the  intensity  of  the 
reconstructed  distribution  at  30®  pitch  angle.  A  plot  of  the  theoretical  diffusion  coefficient 
which  best  fit  the  deconvolved  pitch  angle  distribution  was  selected.  Fits  were  determined 
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Figure  5-5.  Reconstructed  pitch  angle  distribution.  The  solid  line  is  the 
reconstruction  of  the  pitch  angle  distribution  in  1®  increments  produced  by  the 
deconvolution  algorithm.  The  counts  are  plotted  as  crossed  error  bars  where  the 

vertical  error  bar  represents  ac2»  and  the  horizontal  bar  is  the  instrument's  full 
width  response.  The  uncertainty  sgO  is  0.01. 
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Figure  5-6.  Fits  of  observed  electron  distributions  to  theoretically  diffused 
distributions.  The  heavy  solid  line  is  the  deconvolved  distribution  and  the  dashed 
line  is  the  best  fit  predicted  angular  distribution  which  would  result  from  moderate 
(a-c)  to  weak  diffusion  (d-f)  from  Lyons  [1973].  Given  for  each  fit  are  the 
observed  anisotropy  R,  the  deconvolved  anisotropy  Rdec.  the  diffusion 

coefficients  D(0°)  and  D(90®)/D(0®). 
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by  eye.  Figure  5-6  is  an  illustration  of  sample  fits  used  to  derived  the  relative  diffusion  for 
various  energy  levels.  The  heavy  solid  line  is  the  deconvolved  distribution  and  the  dashed 
line  is  the  best  fit  predicted  angular  distribution  which  would  result  from  moderate  to  weak 
diffusion.  Cases  5-6a-c  represent  fits  corresponding  to  moderate  diffusion  while  cases  5- 
6d-f  are  consistent  with  weak  diffusion.  The  original  and  deconvolved  anisotropies  are 
also  noted  in  Figure  5-6. 

Figure  5-7  illustrates  the  various  estimates  of  diffusion  rates  obtained  for  19.4keV 
electrons.  All  the  19.4keV  mean  bin  ratios  in  Figure  5-7a  equated  to  moderate  diffusion, 
save  one.  The  exception  was  the  most  anisotropic  average  ratio  (R=3.20)  located  in  the 
post-midnight  sector  at  6^  L  <7.  This  ratio  was  consistent  with  weak  diffusion. 
Observations  representing  the  mean  anisotropy  in  MLT  (integrated  over  all  L)  and  L 
(integrated  over  all  MLT)  values  are  shown  Figure  5-7b  and  5-7c,  respectively.  Figure  5- 
7d  is  an  example  of  a  nearly  isotropic  distribution.  An  opposite  case,  strong  anisotropy,  is 
shown  in  5-7e  for  the  same  bin.  The  nearly  isotropic  case  corresponds  to  moderate 
diffusion  while  the  very  anisotropic  case  corresponds  to  weak  diffusion.  As  expected 
stronger  diffusion  would  be  required  to  isotropize  the  distribution  in  5-7d. 

A  second  data  set  for  4.52keV  particles  was  subsequently  examined.  An  energy 
level  below  lOkeV  was  suggested  as  being  more  representative  of  the  diffuse  aurora  and 
closer  aligned  to  the  work  of  Lyons  [1974]  [Lyons,  personal  communications,  1990].  To 
get  convergence  a  slight  change  had  to  be  made  to  the  deconvolution  algorithm.  The 
uncertainty  used  in  the  deconvolution  algorithm  is 

^decon  =  ^A2®  ^8®  (5-9) 

where  fx  is  the  electron  count  and  sgO  is  a  constant  The  sgO  factor  widens  the  uncertainty 
to  achieve  convergences  in  cases  with  large  temporal  variations.  The  default  value  of  sgO 
is  0.01.  The  uncertainty  sgO  was  raised  to  0.05  for  several  of  the  4.52keV  electron  angular 
distributions.  The  results  were  six  viable  deconvolved  pitch  angle  distributions  with  weak 
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to  moderate  diffusion  (see  Fig.  5-8).  The  reconstructed  anisotropies  decreased  an  average 
of  25  ±  9%  from  those  obtained  using  only  the  ratios. 

A  final  estimate  of  the  diffusion  rate  was  obtained  using  1.09keV  electron  data. 
This  set  of  data  did  not  meet  the  same  strict  requirements  of  the  previous  ones.  It  was  a 
target  of  opportunity,  having  been  extracted  from  the  same  data  set  used  for  the  4.52keV 
electron  study.  The  anisotropy  of  these  samples  ranged  from  within  1  to  54%  of  the  mean 
anisotropy  of  their  respective  bins  (see  Fig.  5-9).  The  presence  of  an  anti-loss  cone  caused 

the  rejection  of  several  1.09keV  angular  distributions.  The  deconvolved  distributions  had 
an  average  decrease  in  anisotropy  of  21±10%  relative  to  that  obtained  from  the  ratios. 

An  inherent  uncertainty  in  the  manner  in  which  the  raw  counts  were  tabulated  is 
present  in  the  results  shown  in  Figures  5-7,8, 9.  The  pitch  angle  assigned  to  each  second 
of  raw  data  represents  the  center  of  a  6°  angular  bin  (the  angular  distance  covered  by  the 
sensor  in  Isec.).  The  sequential  order  in  which  the  energy  channels  are  observed  plus 
other  instrumental  factors  determines  the  amount  of  pitch  angle  error.  The  greatest  pitch 
angle  errors  would  be  for  the  first  and  last  channel  sampled  within  that  second.  This  meant 
there  was  a  2.8°  error  for  the  raw  l9.4k6V  data  points,  ~0.6°  error  for  the  4.52keV  data 
and  -1.1°  error  for  the  l.()9keV  data.  This  systematic  error  does  not  appear  to  make  a 
significant  difference  in  the  results.  The  estimated  diffusion  rates  for  19.4keV  electrons  is 
of  the  same  magnitude  or  only  slightly  higher  than  that  estimated  for  the  1.09  or  4.52keV 
electrons,  despite  the  -1°  larger  systematic  pitch  angle  offset  for  the  l9.4keV  data. 

5.3  ECH  WAVE  AMPLITUDE 

The  ECH  wave  amplitude  required  to  create  the  reconstructed  distribution  was 
determined  by  adapting  Lyons  [1974]  normalized  pitch  angle  diffusion  coefficient  equation 
(3-47)  and  solving  for  the  wave  amplitude 
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THEORETICAL  DIFFUSION  OF  4.52KEV  ELECTRONS 
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Figure  5-8.  Estimated  diffusion  of  4.52keV  electron  distributions  representing  the 
mean  anisotropy  of  the  quiet  time  bins. 
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THEORETICAL  DIFFUSION  OF  1.09KEV  ELECTRONS 
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Figure  5-9.  (a)  Estimated  diffusion  of  quiet  time  l.OPkeV  electrons  whose 
anisotropy  is  given  in  (b). 


Table  5-3.  ECH  wave  amplitudes. 


CASE  LOCAL  TIME 
(MAG) 

L 

VALUE 

EUECTRON  Rbn 
(keV) 

R 

Rdbcon 

OceV) 

<Doa> 

xl0-> 

Doa  ECH  AMP 
xl(h«  (mV/tai) 

1 

1800-M00 

5 

19.4 

1.78 

1.78 

1.40 

1.4 

2.0 

0.1 

17.8 

2 

2200^00 

5 

19.4 

1.95 

1.96 

1.24 

0.9 

2.0 

0.06 

18.9 

3 

02(XW)600 

5 

19.4 

2.82 

2.70 

1.64 

0.4 

5.0 

0.02 

21.4 

4 

1800-2200 

6 

19.4 

1.88 

1.85 

1.35 

OJ 

1.9 

0.005 

32.1 

5 

2200-0200 

6 

19.4 

2.03 

2.04 

1.42 

0.5 

0.2 

0.04 

5.2 

6 

0200-0600 

6 

19.4 

3.20 

3.44 

2.67 

0.8 

0.2 

0.05 

5.9 

7 

1800-2200 

7 

19.4 

1.86 

1.85 

1.35 

0.6 

1.8 

0.04 

16.2 

8 

2200-0200 

7 

19.4 

2.05 

2.08 

1.79 

1.8 

1.5 

0.1 

16.1 

9 

2200-0200 

6 

19.4 

2.03 

2.70 

1.75 

0.8 

0.2 

0.05 

5.5 

10 

2200-0200 

6 

19.4 

2.03 

1.17 

0.95 

1.6 

2.0 

0.1 

18.5 

11 

2200-0200 

5 

4.52 

2.51 

2.69 

2.30 

0.7 

0.1 

0.2 

2.2 

12 

0200-0600 

5 

4.52 

2.97 

2.69 

1.85 

1.1 

0.1 

0.6 

1.6 

13 

2200-0200 

6 

4.52 

2.26 

2.54 

2.21 

1.3 

0.1 

0.9 

1.3 

14 

0200-0600 

6 

4.52 

2.93 

2.80 

1.85 

0.7 

1.1 

0.2 

6.1 

15 

1800-2200 

7 

4.52 

2.49 

2.70 

2.21 

1.0 

1.0 

0.6 

4.2 

16 

2200-0200 

7 

4.52 

2.13 

2.21 

1.46 

1.1 

0.9 

0.6 

4.2 

17 

0200-0600 

6 

4.52 

2.97 

1.22 

1.04 

1.0 

1.3 

0.4 

5.5 

18 

2200-0200 

5 

1.09 

2.48 

2.30 

1.69 

16 

2.0 

0.48 

19 

020(L0600 

5 

1.09 

2.25 

1.85 

1.48 

1.0 

0.7 

10.0 

0.88 

20 

2200-0200 

6 

1.09 

3.80 

25A 

1.65 

1.3 

0.06 

10.0 

0.28 

21 

1800-2200 

7 

1.09 

2.24 

2.21 

1.77 

0.5 

0.85 

22 

2200-0200 

7 

1.09 

1.69 

1.46 

1.19 

1.1 

0.4 

0.72 

23 

220(H)200 

6 

1.09 

2.19 

1.00 

1.04 

6.2 

9.0 

0.81 

R  -  JA 

Rbin  -  bin  average  R 

Rdecon  *  deconvolved  R 

Kth  -  ihennal  energy 

<Dca>  -  diffusion  cocfTicient 

Daa  -  normalized  diffusion  coefTicienl 

ECH  AMP  -  ECH  wave  amplitude 
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lEwlo^  =  <Daa> - ^ - 

®aa(tto>K/Kth) 


(5-10) 


recall  lE^lo  is  the  equatorial  ECH  wave  amplitude,  <Daa>  is  the  bounce  averaged  diffusion 
coefficient,  K  is  the  particle  energy,  K*  is  the  thermal  energy,  and  iZ\xa(«o>K/Kth)  is  the 
normalized  pitch  angle  diffusion  coefficient 

The  values  for  ©aa  were  available  from  the  calculations  of  Lyons  [1974].  The 
thermal  energy  K*  (actually,  the  average  energy)  was  derived  using  an  algorithm  which 
integrates  the  entire  21  channel  energy  distribution  in  a  three  dimensional  model.  This 

requires  an  assumption  of  symmetry  about  the  parallel  axis  [Croley,  personal 
communications,  1990].  The  thermal  energies  used  in  (5-10)  were  an  average  of  Kthj^  and 

Kth|| .  In  Lyons  [1974],  and  Kthu  were  considered  equal.  The  case  to  case  variation 

of  Kth  in  the  19.4keV  sample  cases  were  large,  Ku,  =  0.83  ±  0.52.  Much  smaller  case  to 
case  variations  of  thermal  energy  were  noted  in  the  sample  cases  used  for  4.52  and 
1.09keV;  K*  =  0.99  ±  0.2  and  1.1  ±  0.2,  respectively. 

The  bounce  averaged  diffusion  coefficient  <Dao>  was  derived  from  the  results  of  the 
fit  performed  to  determine  the  strength  of  the  diffusion.  The  minimum  lifetime  (3-43) 
was  factored  out  of  the  diffusion  coefficient  obtained  from  Figure  3-2  so  that  <Daa>  = 
D(0°). 

The  resultant  ECH  wave  amplitudes  represent  the  minimum  wave  emissions 
required  by  the  theoretical  results  of  Lyons  [1974]  to  produce  a  distribution  similar  to  the 
observed.  For  the  19.4keV  electrons  in  Table  5-3,  the  ECH  wave  amplitudes  ranged  from 
5.2  to  32.1mV/m.  Two  cases  with  anisotropies  much  higher  and  much  lower  than  the  bin 
average  (Rbin=2.03)  in  the  midnight  sector  at  6<  L  <7  were  also  determined.  For  a  very 
anisotropic  distribution  (R=2.70),  case  9  on  Table  5-3,  the  required  wave  amplitude  was 
5.51  mV/m.  Case  10  is  a  nearly  isotropic  distribution  (R=1.17)  and  would  require  an 
amplitude  of  18.5mV/tTi  This  change  in  wave  amplitude  reflects  the  greater  power  case  10 
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requires  to  isotropize  the  distribution. 

Six  cases  of  average  anisotropies  were  available  at  the  4.52keV  energy  level.  The 
estimated  wave  amplitudes  corresponding  to  these  4.52keV  angular  distribution  are  cases 
1 1  through  16  on  Table  5-4.  They  ranged  from  1.3  to  6.1mV/m.  A  nearly  isotropic  case 
(R=1.22)  in  case  17  would  be  consistent  with  a  5.5  mV/m  wave  amplitude.  This  is  a 
relatively  high  amplitude  and  is  indicative  of  the  greater  power  required  to  fill  the  loss  cone. 
Wave  amplitudes  of  0.28  to  0.88mV/m  are  needed  to  reproduce  the  anisotropies  of  the 
1.09keV  electron  distributions.  The  wave  amplitudes  consistent  with  the  anisotropy 
observed  in  each  local  time/L  bin  are  cases  18  through  22  on  Table  5-4.  The  isotropic 
sample  (R=1.00)  in  case  23,  as  noted  before  in  similar  cases  above,  would  require  an 
relative  high  amplitude  of  0.81mVAn  to  create  the  observed  distribution. 


67 


CHAPTER  6 

CONCLUSIONS 


I  have  investigated  the  assumption  that  electrons  in  the  near  earth  plasma  sheet  are 
on  strong  diffusion  during  magnetically  quiescent  times.  This  assumption  has  been 
invoked  in  numerous  studies  of  magnetospheric  coupling  with  the  diffuse  aurora  [Fontaine 
et  al.,  1985,  Fontaine  and  Blanc,  1983;  Kennel  and  Ashour-Abdalla,  1982;  Kennel,  1969; 
among  others].  The  assumption  of  strong  diffusion  is  supported  by  low  altitude 
observations  of  isotropy  on  precipitating  electrons  causing  the  diffuse  aurora 
[Winningham,  1975;  McDairmid  and  Budzinski,  1968;  Chase,  1968]  and  isotropy  of  the 
plasma  sheet  electron  distributions  during  geomagnetically  active  periods  [Baker  et  al., 
1981].  In  this  study,  sample  pitch  angle  distributions  were  selected  to  represent  average 
electron  angular  distributions  for  the  quiet  time  (Kp<3)  plasma  sheet.  The  instrument 
response  was  removed  from  the  sample  distributions  by  means  of  a  deconvolution 
algorithm. 

A  total  of  nineteen  cases  representing  average  pitch  angle  anisotropy  were  examined 
at  three  energies  (1.09,  4.52,  and  19.4keV).  None  of  the  cases  displayed  angular 
distributions  consistent  with  strong  pitch  angle  diffusion.  Moderate  to  weak  diffusion  was 
found  exclusively.  Examination  of  a  case  of  nearly  isotropic  distribution  at  19.4  and  one  at 
1 .09KeV  also  were  consistent  with  only  moderate  diffusion.  I  thus  conclude  that  electrons 
in  the  equatorial  near  earth  plasma  sheet  are  usually  not  under  the  influence  of  strong 
diffusion  when  Kp^. 

The  assumption  that  strong  diffusion  is  required  to  sustain  the  diffuse  aurora  is  then 
unsupported.  This  lack  of  strong  diffusion  is  inconsistent  with  the  previously  observed 
isotropic  pitch  angle  distributions  in  the  diffuse  aurora.  One  possible  explanation  is  that  a 
secondary  process  is  acting  to  isotropize  the  precipitating  electrons  at  low  to  middle 
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altitudes.  Koons  et  al.  [1972]  suggested  scattering  by  extremely  low  frequency  waves 
observed  at  altitudes  of  20(X)-4500km  might  be  responsible.  Another  possibility  is  that  the 
pitch  angle  isotropy  observed  at  low  altitudes  in  the  diffuse  auroral  regions  is  not  the 
predominate  distribution  during  geomagnetically  quiet  times.  A  review  of  low  altitude 
electron  data  limited  to  magnetic  quiescence  is  suggested  as  a  follow-up  study. 

To  determine  if  weak  to  moderate  diftiision  evidenced  by  the  results  is  adequate  to 
create  and  sustain  the  diffuse  aurora,  the  flux  in  the  atmospheric  loss  cone  must  be 
checked.  This  would  require  a  comparison  of  equatorial  electron  flux  in  the  loss  cone  to 
the  precipitating  flux  in  the  diffuse  aurora.  Simultaneous  observations  from  high  and  low 
altitude  satellites  near  the  same  magnetic  field  line  would  be  required.  To  derive  the  flux  in 
the  loss  cone,  either  the  resolution  of  particle  detectors  must  be  improved  or  observations 
must  undergo  an  extensive  deconvolution  process  such  as  was  done  in  this  study.  Also, 
great  care  must  be  taken  with  the  magnetic  mapping. 

The  second  part  of  this  investigation  dealt  with  the  minimum  ECH  wave  amplitudes 
necessary  to  create  the  observed  pitch  angle  distributions.  All  other  possible  contributors 
to  the  diffusion  process  were  excluded  from  this  study.  The  resulting  wave  amplitudes 
ranged  from  a  quiet  time  high  of  32.1mV/m  for  the  19.4keV  electrons  to  a  low  of  .28mV/m 
for  1.09keV  electrons.  These  ECH  wave  amplitudes  are  one  to  three  orders  of  magnitude 
greater  than  those  observed  by  Roeder  and  Koons  [1989]  and  Belmont  [1983]  (see  Table 
6-1).  These  high  amplitudes  are  in  agreement  with  the  results  of  Schumaker  et  al.  [1989] 
who  also  derived  higher  than  observed  amplitudes  from  their  set  of  equatorial  electron 
distributions.  Recall,  the  distributions  reported  here  were  found  to  be  undergoing  only 
weak  to  moderate  diffusion.  If  strong  diffusion  were  in  effect  the  derived  amplitudes 
would  have  been  higher  by  a  factor  of  3  to  10. 

Thus,  the  conclusion  that  particle-ECH  wave  interactions  are  the  primary 
mechanism  diffusing  electrons  into  the  loss  cone  at  during  quiet  times  is  rejected.  The 
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Tabic  6-1.  Occurrence  of  ECH  wave  emissions  [From  Rocder  and  Koons,  1989J. 


question  of  what  wave-particle  interactions  or  other  process  is  responsible  for  the 
continuous  diffuse  aurora  remains  unanswered.  For  lack  of  any  other  known 
semicontinuous  source  of  waves  in  the  equatorial  plasma  sheet,  Roeder  and  Koons  [1989] 
suggested  the  possibility  that  wave-pardcle  interactions  at  low  or  middle  altitudes  on  auroral 
field  lines  might  be  responsible  for  the  near  pitch  angle  isotropy  observed  there  in 
conjunction  with  the  diffuse  aurora. 

I  hope  the  observations  presented  here  will  aid  in  the  continuing  search  for  the 
elusive  mechanism  coupling  the  ionospheric  diffuse  aurora  to  the  magnetosphere. 
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